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Interdisciplinary Approach to Propulsion System Simulation 

The "computational test-cell" will enable the Incorporation of new 
methodologies, such as concurrent engineering and probabilistic methods, into 
the propulsion design process. This will provide the capability to conduct 
credible, Interdisciplinary analyses of new propulsion concepts and designs. 

Probabilistic methods can be used as the basis for reliability-based 
design. Recently methods have been devised that provide the capability of 
simulating the performance of propulsion systems at several levels of resolu- 
tion. These methods make It possible to quantify uncertainty and to establish 
confidence bounds for the calculated values. 

The Introduction of reliability-based design methodology along with 
probabilistic analyses will provide a tool to reduce the design space for new 
systems and to reduce our dependence on hardware testing for proof-of-concept 
and system Integration demonstrations. The resulting simulations will reduce 
the need for testing and Identify potential operational problems early In the 
design process. 

This capability will make It possible to compute the expected perfor- 
mance, stability, reliability, and life of propulsion components, subsystems, 
and systems at design and off-design conditions, to bring life cycle cost 
trade-offs early Into the design process and to determine optimum designs to 
satisfy specified mission requirements. 

Physics Modeling 

The analysis of propulsion phenomena Involves a combination of 
disciplines Including fluid mechanics, thermal sciences, structural mechanics, 
material sciences, combustion theory and controls theory. The degree of 
resolution within an analysis Is determined by the magnitude of local effects, 
the extent of their region of Influence, and the dynamic time scales of the 
appropriate physics relative to the dynamic scale of the system phenomenon 
being analyzed. Often the limiting factor will be the available computer 
power (speed and memory capacity). The analyst must determine which terms In 
the governing equations to retain and which to Ignore so as to achieve the 
maximum level of fidelity within the computational constraints. 

The coupling of the disciplines and component codes Involves the 
subdivision of a complete system, e.g., an aircraft engine. Into a series of 
subsystems, e.g., Inlet, compressor, combustor, turbine, and nozzle. It Is 
convenient to define a hierarchy of multi -disciplinary simulation modules for 
each subsystem ranking from relatively simple time and space "averaged" 
analysis methods (level I) to complex three-dimensional, time-accurate 
analysis methods (level V). 

Level I: Engine system performance model. This model Is basically a 
thermodynamic model which calculates the system efficiency based upon 
engine configuration and component efficiencies. It allows rapid 
evaluation of various engine concepts. 

Level II: Engine system dynamics and controls model. This model Is 

basically a one dimensional flow path model, with simplified structural 
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elements, controls, and other disciplines. It uses component 
performance information, design geometry Information, and dynamic 
information In order to calculate engine thrust and weight as well as 
system transient response In order to analyze operability problems and 
devise control strategies to handle them. 

Level III: Space and/or time-averaged engine system model. This model 
Is basically a two-dimensional (l.e., ax 1 symmetric) fluid model. It 
utilizes axisymmetrlc, coupled discipline models In an engine system 
environment In order to relate component boundary conditions (primarily 
Input/output conditions) to overall system boundary conditions in order 
to simulate component Interactions. This Is also the basic level about 
which the "zooming* process Is constructed. It will be a transient 
model and address all problems from level II but, in addition, provide 
more detailed geometry information. 

Level IV: Space and/or time averaged sub-system (or component) models. 
These models are basically three dimensional. They are multi-discipline 
models which are coupled in ways which are compatible with the physics 
of the component model, but are still averaged over smaller time and 
space scales. These models must also be post-processed in order to 
connect with the level three engine system model in the "zooming" 
process. 

Level V: Three-dimensional, time-accurate component models. This level 
of simulation basically consists of a fully three-dimensional, time- 
accurate simulation of all physical processes on a component-by- 
component basis. This Is the most complete level of physical 
approximation. 

Dlcslpllne Coupling 

Propulsion phenomena are inherently multi-disciplinary (l.e., the true 
system response Is the coupled effect of all the participating disciplines and 
the aggregate of the system components' responses and interactions). Present 
analyses (and experiments) tend to focus on single-discipline aspects of the 
phenomena within a local region (e.g., a single component). Using suitable 
approximations, these analyses are sometimes extended to a propulsion 
subsystem or, In rare cases, the complete propulsion system. 

Recent advances In computational fluid mechanics, computational 
structural mechanics, computational materials science, computational controls, 
and computer science and technology make it feasible to consider the 
development of a "computational test-cell" for propulsion that would allow for 
comprehensive simulation and analysis of entire propulsion concepts and 
designs before committing to hardware. 

For computational simulation to be credible, it must include efficient 
multi-disciplinary coupling. In the case of multi-disciplinary simulation of 
dynamic phenomena, the time scales associated with various aspects of the 
phenomena have to be considered. In an engine Interacting phenomena, such as 
surge, stall, flutter, component and system dynamics, low and high cycle 
fatigue, and takeoff and landing operations occur within widely varying time 
intervals. The computational procedures and the "clock cycle" of a multl- 
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disciplinary simulation have to accommodate these vast differences In time 
scales. The simulation clock cycle has to be consistent with the available 
computational power and. In the case of animated graphic representation, the 
perception rate of the human visual capability. 

Implementing coupling In the required numerical simulation, analysis, 
and optimization Is a tremendous challenge because of the potentially very 
large number of Interrelated variables and the very large number of Iterations 
that can result from general-purpose algorithms. A hierarchical approach that 
can reduce the dimensionality of the system description while still retaining 
the essential system behavior is needed. There are a variety of techniques 
that can be used for coupling discipline variables for propulsion components, 
subsystems, and systems. These Include sequential Iteration between 
disciplines, specially-derived system matrices, and coupling at the 
fundamental equation level. In NPSS, all three methods will be applied to the 
filtered Navler-Stokes equations and the progressively substructured 
structural mechanics formulations. Relationships (l.e., sensitivities) will 
^? r use optimization algorithms that are streamlined for the 
multi-disciplinary, multi -component application. 

The coupling across disciplines In a concurrent multi-disciplinary 
formulation can be represented by coupling relations. The coefficients 
(elements) In these relations define the coupling of a specific variable from 
one discipline with respective variables from Interacting disciplines (Fig. 

6). Perturbation of the variables in the coupling relations provides a 
measure of the sensitivity of the Interacting disciplines to this 
perturbation. A priori description of this sensitivity relationship enhances 
the computational simulation In several respects: 1) scoping the degree of 
coupling, 2) Identifying the Interacting disciplines, 3) resolving time/space 
scales, 4) selecting time/space scale for loosely coupled Interacting 
discipline Intervention during the solution processes, 5) deciding on a 
solution strategy, and 6) imposing convergence criteria. 

Analysis Fidelity 

*u A ^ e, "P tin 9 to resolve all of the length and time scales that are present 
In the fluids and structures of the engine Is impractical, even on high 
performance computers. Therefore, a rational approach for identifying and 
resolving the critical scales Is needed. Approaches that have been shown to 
be effective for single component analyses will be extended to the simulation 
of coupled components and entire engine systems. Approaches will be developed 
to allow selected components to be resolved to a greater level of detail than 
others, Utilization of the zooming approach will allow the Interconnection of 
a series of multi-discipline simulations In which a single or small number of 
modules are simulated with very accurate methods, perhaps level IV or V, while 
the remainder of the subsystems are Implemented with simple methods, perhaps 
level II or III. This focusing or "zooming" In on a particular component will 
anow for a more thorough analysis of that subsystem In a complete multi - 
discipline system format without having to completely simulate the entire 
system at the same detailed level. 

, * For example, studies of compression system stability will require a 
detailed treatment of the compression system to be coupled to lower-resolution 
treatments of the fan, combustor, turbine, and nozzle with the appropriate 
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boundary conditions to represent the inter-component interactions (Fig. 5). 
This "zooming" capability will permit the analyst to capture relevant physical 
processes throughout the engine In a computational ly-tractable manner and will 
allow the analysis to be used on a routine basis for design assessment and 
optimization. Thus, this approach will be much more cost effective and should 
provide an attractive approach for overall system performance optimization. 

The actual interface algorithms used in this zooming approach will range from 
the direct coupling approach described above to one Involving the interface of 
time- and space-averaged parameters. With this approach, special emphasis can 
be placed on the effects of Interface sensitivities between two subsystems in 
an entire system. 

Computational Simulation 

Digital simulation of aerospace propulsion system behavior has been In 
existence for many years. The earliest simulations were developed In the 
seventies. The performance and reliability of engine systems depend on the 
dynamic Interaction of their subsystems which. In turn, depend on the dynamic 
interaction of their respective components. Interaction phenomena of 
importance Include flutter, rotor Instability, fatigue, flow separation, 
nonuniform combustion, blade containment, and noise suppression. The determi- 
nation of aero-thermodynamic system performance has traditionally relied on 
prototype tests while structural reliability has been calculated from field 
data. This experience has been used to develop simulation techniques that 
employ varying degrees of approximation to model and compute the aero-thermo- 
dynamic performance and structural reliability of new designs. In general, 
these simulations can be divided into two classes, depending upon the time 
dependance. 

Steady state simulations are normally used by design engineers In order 
to assess design trade-offs. Here the emphasis Is on ease-of-use by the 
designer and, In particular, allowing the designer to Include "company lore" 
or company expertise In the design. Depending upon the use of the design 
system (I.e., whether It Is for conceptual design, preliminary design, 
detailed design, or final design) there will be more or less fidelity included 
in the simulation. Steady-state simulations are used for design points 
analysis, with allowances for off-design performance. In the latter design 
stages, steady-state simulations can be used to develop control system 
schedules and to provide estimates of engine system life. 

Dynamic simulations are used after the engine Is designed In order to 
develop control laws/logic and to determine the limits of stable engine 
operation. Obviously, If the simulation calculations can be speeded up, more 
detail (I.e., spatial and temporal resolution) can be Included In the 
simulation model. During control system hardware and software Implementation, 
there is a need for a real-time engine simulation that can be operated with 
the control system In a "closed-loop" fashion. 

Dynamic simulations are also used to study cases when the engine behaves 
differently In the field from what was envisioned In the design phase, or as 
uncovered In the testing of the engine before It was installed Into service. 
These simulations can be particularly valuable when ground-based experimental 
facilities are not available to simulate the In-flight conditions under which 
the unusual behavior was observed. Obviously the more accurate the 
simulations, the more their value. 
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The development of an engine simulation capability will begin with 
existing Level II dynamic engine system models of aerothermal and structural 
behavior. Level IV aerothermal and structural simulations will be used to 
generate the required component parameters and maps for the Level II engine 
models. Then, methods for Improving the parametric representation of the 
components will be Investigated so that the significant phenomenon observed 
from detailed analyses can be represented In the engine model. 

The Initial simulations Involving Level II aero and structures codes 
will Investigate the thermal lag between changes In the engine operating 
conditions and the heat transfer effects on the structure. Thermal strains 
resulting from the changes In the temperature of the structure affect the 
secondary cooling flow passages and tip clearance flow In the components. 
These effects must be accounted for in the aero codes and will result in a 
change In the computed engine operating conditions. 


The Level IV aerodynamic simulation model that will serve as the basis 
for the Integrated propulsion system model will be the Adamczyk average- 
passage formulation which consists of the filtered forms of the Navler-Stokes 
and energy equations. This model was designed to resolve only those temporal 
and spatial scales that have a direct impact on the relevant physical 
processes. The effects of the unresolved scales, which appear as body forces 
and energy sources In the equations, are estimated through semi -empirical 
relations, based on experiments or high -resolution numerical simulations. The 
results from the lower- resolution analysis appear as boundary conditions for 
the high-resolution model. Initially, this model will be applied to the study 
of a compression system and Its performance, stability, blade vibration, and 
noise generation. Since the methodology applies to the fundamental fluid flow 
equations, It will then be extended to the other propulsion components. 


The structures modeling will be aimed at developing a comparable 
computational capability that will provide a means to traverse multiple scales 
of spatial resolution with a minimum number of variables at each level In 
this way, an analysis can proceed from a blade to a rotor sector to a rotor to 
an engine core to the complete engine. The resulting system model will have a 
minimum number of degrees of freedom consistent with the objectives of the 
analysis which will minimize the computational requirements. This methodology 
will be applicable to the solution of linear and Incremental nonlinear 
analysis problems. This capability will be achieved through the formulation 
and Implementation of a progressive substructuring ("telescoping super- 
elements") technique within the mixed-iterative finite element method 
framework and associated MINUTES computer code and within the boundary element 
framework and associated BEST3D computer code. 


Computing Platform Portability 


It Is the Intention of the NPSS to take advantage of existing codes to 
the extent possible, while at the same time maintaining the flexibility to 
utilize emerging massively parallel computing hardware platforms. The 
architecture envisioned ustlllzes shared memory programming paradigm and 
standard software tools and programming language extensions. The programing 
will be Independent of hardware architecture. 


Within the computing system, the nature of the coupling between the 
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computing system components (processor, memory, communication) will depend 
upon the engine system component codes and single discipline codes required to 
compute the desired engine attributes. Therefore, the selection/development 
of appropriate processor I/O software, compilers, networking protocols will be 
accomplished In conjunction with the development of engine system and 
discipline (l.e. application) codes. 

The NPSS technology project can develop and demonstrate many key, 
enabling technologies for aerospace propulsion systems design, analysis, and 
optimization. However, to be successful, several things must take place: 

Effective Inter-disciplinary teams must be established to define, 
advocate, and Implement technical solutions. 

Coordination and a balancing of efforts among the Inter-disciplinary 
engine system activities (physics, algorithms, models, codes) and the 
Inter-disciplinary computer system activities (architectures, software 
tools,..) must be maintained so as not to push either activity ahead of 
the other. 

This suggests a strong requirement for effective project management to 
ensure that the available funding and skilled staff are effectively used to 
address the needs. 
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NPSS Goal 

• Reduce life-cycle costs by advancing system analysis 
capability through high-fidelity computational 
simulations by 

- Higher level of concurrent engineering 

- Rapid evaluation of effects of new and novel 
concepts on system performance 

- Early risk assessment 

- Early operability studies 

- Rapid evaluation of field problems 

- Assessment of performance degradation 


NPSS LEVELS OF MODELING FIDELITY 


LEVEL 1 : ONE DIMENSIONAL STEADY-STATE 


Oita Dimensional Engine System Analysis 
Performance Maps 
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LEVEL 2: ONE DIMENSIONAL TRANSIENT 


On* Dimensional Transient Engine System Analysis 

Psrtormance Maps Mixing Volumes 
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dt 


LEVEL 3: AXISYMMETRIC/TWO DIMENSIONAL 
QUASI-STEADY-STATE 


Axisym metric Engine Analysis 


LEVEL 4: THREE DIMENSIONAL 
QUASI-STEADY-STATE 


Three Dimensional Engine Component Analysis 
Steady Flow Sink Source 


Fan 

Analysis 




LEVEL 5: THREE DIMENSIONAL 
TRANSIENT 


Three Dimensional Tranaiant Engine Component Analysis 
Unsteady Flow Sink Sourc# 


Fan 

Analysis 
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Coupled-Multi-Discipline Representation for 
Aerospace Propulsion Systems 



System 

response 

variables 



System definition/characteristics 


and coupling relationships 


A Aero F Fabrication 

T Thermal P Performance 

S Structural C Cost 

M Material 



System 
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and service 
parameters 


NPSS Models 


Acf0 Structures 
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Koen (1985): "...it is the engineering method or design process, rather than the 

artifacts designed, that bind all engineering disciplines together 
and defines the engineer." 


The Engineering (Design) Process 



Deductive Analysis 


Engineering (Designers/Synthesizers) = Engineering Science Technicians (Analysts) 

• Industry Needs Both • 
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